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1. Introduction 

The localization of energy for gravitating systems remains an outstanding problem. 
Unlike all other source and interaction fields, the standard techniques for identifying 
an energy-momentum density for the gravitational field have only yielded various 
expressions which are inherently reference frame dependent; these non-covariant 
expressions are often referred to as pseudotensors. This non-covariant feature can be 
understood as an inevitable consequence of the equivalence principle, which precludes 
the detection of the gravitational field at a point — so one cannot have a point-wise well- 
defined energy-momentum density for gravitating systems (for a good discussion of this 
point see Ch. 20 in [1]). 

The energy-momentum pseudotensor approach has largely been displaced by the 
more modern perspective of quasi-local: energy-momentum is to be associated with a 
closed 2-surface (for a review of the quasi-local idea see [2]). One quasi-local formulation 
is in terms of the Hamiltonian. The Hamiltonian for evolving a (generally finite) 
spacetime region includes a boundary term. Quasi-local quantities are associated with 
this Hamiltonian boundary term. There are many possible quasi-local expressions 
simply because there are many possible boundary terms. They are all physically 
meaningful, for each distinct boundary term is associated with a distinct physical 
boundary condition (which is given by what must be held fixed in the Hamiltonian 
variation). We note that this Hamiltonian quasi-local approach includes all the 
traditional pseudotensors, since each is generated by a superpotential which serves 
as a special type of Hamiltonian boundary term [31 HI |7]. Via the Hamiltonian 
formulation the ambiguities of the traditional pseudotensors (which expression? which 
coordinate system?) are clarified. Hence, from the perspective of this Hamiltonian 
boundary term approach to quasi-local energy-momentum the traditional pseudotensors 
are still of interest. 

For pseudotensor expressions the physical energy-momentum of the gravitational 
field is inextricably bound up with the choice of reference frame. This problem has long 
been recognized. Indeed soon after Einstein had proposed his expression for gravitational 
energy it was noted that with certain choices of coordinates it could give both a nonzero 
energy for Minkowski space and a vanishing energy for the Schwarzschild solution P,|9]. 
Notwithstanding this, in many cases the reference frame ambiguity is not a problem — 
simply because there is an obvious natural choice of reference frame. 

In particular in evaluating the energy-momentum for an asymptotically flat 
gravitating system, the asymptotic Minkowski space provides a natural and 
unambiguous reference for the coordinate system. For this case almost all the classic 
pseudotensors (except for M0ller's 1958 expression [10]) give the standard value for 
the total energy-momentum. Now there is another case where there is a natural and 
unambiguous reference: namely in a small region, where one can use the flat tangent 
space at some interior point to determine a Minkowski coordinate system. Here we shall 
test the classical pseudotensors in this limit. 
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Specifically we will be concerned with the famous pseudotensors due to Einstein 
[IT] . Papapetrou [12], [131 EH] > Landau-Lifshitz [15], Bergmann-Thompson [IB] . Goldberg 
[17] . M0ller [ID], and Weinberg [181 HJ- Although in some interesting cases many of 
the pseudotensors give identical answers (see [IH]), in the small vacuum region limit 
considered here that will not be the case. 

Note that a good energy-momentum expression for gravitating systems should 
satisfy a variety of requirements (see, e.g., [2Ql[2]), including giving the standard values 
for the total quantities for asymptotically flat space and reducing to the material energy 
momentum in the appropriate limit (equivalence principle). No entirely satisfactory 
expression has yet been identified. One of the most restrictive requirements is positivity. 

It is generally accepted that gravitational energy should be positive; indeed positive 
energy proofs have been heralded (e.g., [21] [22], [20]). Positivity is difficult to prove 
in general. One can regard positivity as an important test for quasi-local energy 
expressions. One limit that is not so difficult, and which has not been systematically 
investigated for all the classical pseudotensor expressions, is the small region limit. 
The small region requirements have not yet been applied to many energy-momentum 
expressions. We found that they afford both interesting restrictions and unexpected 
freedom. 

For a small region within matter the equivalence principle requires that the energy- 
momentum expression should be dominated by the material energy-momentum tensor. 
On the other hand the positivity of gravitational energy in a small vacuum region is 
assured if its Taylor series expansion in Riemann normal coordinates is at the second 
order a positive multiple of the Bel-Robinson tensor. 

Some time ago Deser, Franklin and Seminara [23J presented a discussion of 
pseudotensors in a small region, which has been a major sources of inspiration for us. In 
that work the main techniques that we will use here were developed. They examined the 
Taylor expansion of a pseudotensor around a preselected (vacuum) point using Riemann 
normal coordinates (RNC). Considering pseudotensors derived from superpotentials, 
it was noted that the leading order non-vanishing vacuum expansion was at second 
order and would be quadratic in the Riemann tensor. They then identified four basis 
expressions for such terms and related them to the famous Bel-Robinson tensor. 

They argued that it is desired that one get the Bel- Robinson tensor, B a p^ v . 
They found exactly one such expression from a certain linear combination, B a p^ u = 
d 2 ^ u (L a f3 + \E a p), of the Landau-Lifshitz and Einstein pseudotensors, a combination 
which they argued was unique. 

We have reexamined the issue, considering all the aforementioned pseudotensors. 
Here we report in detail on our results, which were first obtained in |24J and briefly 
announced in [25] . We found that no classical holonomic pseudotensor gives the required 
second order vacuum result (surely this contributed to the difficulty in finding a positive 
gravitational energy proof). However, using similar methods to the work just mentioned, 
we found another independent combination of the Bergmann-Thompson, Papapetrou 
and Weinberg pseudotensors, and thus a one parameter set of pseudotensors, with the 
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same desired Bel-Robinson property. (This was overlooked in the earlier work because 
they had not allowed for pseudotensors which require the explicit use of the Minkowski 
metric as a reference). 



2. The classical pseudotensors 

The classical pseudotensors can be obtained by suitably rearranging Einstein's equation: 

Gjxu = kT^,,, (1) 

(where k := 87rG/c 4 ). One specific way is to choose a suitable superpotential U^ x u = 
U^ x \ and define the associated gravitational energy- momentum density pseudotensor 
by 

2Kt^:=d x U^ v -2\g\^G^. (2) 

Einstein's equation now takes the form 

2kT\ := 2K(\g^T\ + t\) = d x U x » v . (3) 

Because of the antisymmetry of the superpotential the total energy-momentum density 
complex is automatically conserved: d^T^y = 0. 

There are some variations on the idea, the classical pseudotensorial total energy- 
momentum density complexes all follow from associated superpotentials according to 
one of the patterns 

2kT^ u = 8 x U^ u , 2kT^ = x U^, 2kT^ = d aP H a ^\ (4) 

where the superpotentials have certain symmetries which automatically guarantee 
conservation: specifically U^y = U^ u , U^ Xu = U^ u , while H a ^ y has the symmetries 
of the Riemann tensor. In particular the Einstein total energy-momentum density 
follows from the Freud superpotential [26] 

u$\ := -M (5) 

while the Bergmann-Thompson [16], Landau-Lifshitz [15], Papapetrou [12], Weinberg 
[T8] and M0ller [10] expressions can be obtained from the respective superpotentials 



TT\1XV 

<_7 BT 



g u5 ug x s, (6) 

\g\^U&\ equivalents Hff v := MWV™, (7) 

OSrtblVn (8) 
hT v ■= sz&g\k ab (-9 mc g nd + \g mn g cd )g cd , (9) 

U£ v ■= -\g\*g^T a p v 5£ = \gfig^g xs {d mu -d mi/ ). (10) 

(Note that all indicies in this work refer to spacetime and range from to 3, otherwise 
our conventions follow pp.) 

These expressions are all non-covariant. As they depend on the coordinates in 
a non-tensorial way, they can at best be expected to give sensible energy-momentum 
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values only in certain coordinates — which are in some suitable sense nearly Minkowski 
coordinates. Given that we have such coordinates we also have an underlying Minkowski 
space reference structure. Some of the superpotentials explicitly include this reference 
metric, which here has the Minkowski values = diag(— 1, +1, +1, +1). Our basic 
philosophy is that energy-momentum is properly a covector, and hence properly the 
energy-momentum density should be a weight one density with the index positions #V 
For all of the classical pseudotensors this can be achieved by introducing suitable factors 
of the Minkowski metric and its determinant. According to this philosophy only the 
Einstein and M0ller expressions have proper expressions that do not explicitly need the 
Minkowski metric associated with the chosen coordinates. 

It should be noted that, while in general identifying physically meaningful 
Minkowski coordinates is problematical, in the small region case of interest to us here 
at any chosen point there is a natural local Minkowski structure, as we discuss in the 
next section. 



3. Some Technical background 

3.1. Riemann normal coordinates (RNC) 

As Riemann first argued (see, e.g., Spivak [27]). at any preselected point one can choose 
coordinates such that at the point x^ = 0, the metric coefficients have the standard flat 
values, the first derivatives of the metric vanish, and the second derivatives have the 
minimum number (20 for n — 4) of independent values. Specifically 

9a/3\o — 9a.pi d^g a p\o = 0, 3dp V g a p\o = —(R al j,j3 U + Ravpfjjlo, (11) 

where R a is the Riemannian curvature tensor, and, in our case, g a p = r\ a p = 
diag(— 1, +1, +1, +1) is the Minkowski spacetime metric. The corresponding Levi-Civita 
connection values are 

r Q /3 7 |o — 0, 3<9 M r a /3 j,|o = —(R a f3vn + R a vf3[i)\o- (12) 



3.2. Quadratic curvature basis 

It turns out that when expanded in RNC the lowest non-vanishing vacuum energy- 
momentum expressions are of the second order and are quadratic in the curvature 
tensor: t^ v ~ {R....R .' .)^ U ijX l xK An investigation [23j of all such possible terms (taking 
into account the Weyl = vacuum Riemann tensor symmetries) shows that they can be 
written in terms of 

fiauf3 ■ Ra^ibaR v (3 = Quf}^icx = Qafif3ui (l 3) 

where all the symmetries have been indicated. The cited work defined the three basis 
combinations 
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along with the trace tensor 

Tfiva/3 = — ~g p U Q° ' a /3a ■ (15) 

One could write all our vacuum expansions to second order as linear combinations of 
X, Y, Z, T. However it is more suitable for physical purposes to use the Bel-Robinson 
tensor. 

3.3. The Bel-Robinson and two other tensors 
The Bel-Robinson tensor 

B^uafB RppcraR^ u f3 RppaaR^ u {3 pyRaparRff (1^) 

has many well known remarkable properties, sec e.g., [23, 28J. For our considerations we 
are interested in it only in the vacuum, where the Riemann tensor reduces to the Weyl 
tensor. In this case the Bel-Robinson tensor is completely symmetric and traceless, and 
the last term admits an alternate form using 

4:R ap(TT Rp f " TT = gapRnX^&R^ 5 ■ (17) 

In vacuum the Bel-Robinson tensor B and two other convenient tensors S and K are 
given by 

B a /3pu '■= RaXfiaR/3* ' v° + R-aXua Rf3* ~ ~ 9a/39 p.vRnX^R^' 76 , (18) 

o 

S a f3pv '■= RapXaR/3u Xa + RauXaRpp,^ + -^gaf3g ^uRkX-isR^^ , (19) 

3 

Kaf}p,v '■= RaXf3aRp, X v U + RaXfiaRu^ p ~ -QclPQp.vRk.X'i&'R*'^ '• (20) 

O 

In term of the aforementioned basis for quadratic terms 

B = Z + 3T, S = -IX + 2Z - 6T, K = Y + 9T. (21) 

(Here and below we have suppressed some obvious indicies.) 

Our leading order non-vanishing vacuum expressions will be given as linear 
combinations of B, S, and K. As can be directly verified, these three combinations 
(by virtue of (ITT)) ) satisfy the divergence free condition, 

dp^xHi^) = 2x j t ij ai = 0, (22) 

and all such tensors are some linear combination of these three. While the tensor S has 
been known for a long time, we here draw attention to the tensor K which also enjoys 
this divergence free property. 

4. The small region limit 



Here the small region values for the classical pseudotensors are presented. 
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4-1. Einstein 

The Einstein total energy-momentum complex can be obtained from the Freud 
superpotential (jSJ): 

2kT£ v := d x Uf v := d x {-\g\^g^Y a Pl 8±) (23) 

= - l<?l V'C^Vt - r a 6X r s ,)8± - d x (\ g \h Pa W a ^± (24) 
= 2\ g \*G» v - |<^[r^ A r- 7 - r<^ A rvK^. (25) 

Using the Einstein field equation (jTJ this relation takes the form 

T£ u = \g\ l iT» u + t£ u , (26) 
where the Einstein energy-momentum pseudotensor density is 

C := -M^M^r^ - r^ A r^ 7 ]^. (27) 

Expanding in normal coordinates, to zeroth order the pseudotensor vanishes, so (1261) 
reduces to T^ u = |g|^T M „, which is the matter interior limit expected from the 
equivalence principle. In vacuum 7g = t^, and the first non-vanishing contribution 
appears in the second order (note: all trace T terms are proportional to Ricci and hence 
vanish in vacuum); after some computation we find in vacuum 

2 K T£ V = 2Kt% = |^(-2r^ Q rv + ^ a r%) ~ ^y^( 45 ~ S )*A> ( 28 ) 

a result known for some time [TJ [221 [221 [2H] • 

4-2. Bergmann- Thompson, Landau- Lifshitz and Goldberg 

The total energy-momentum complex proposed by Bergmann and Thompson [16J can 
be obtained from a transvected version of the Freud superpotential ([6]): 

2*1% := dxisT'U^s) = 2KT£ s g vS - [glh^ ^& x g" s . (29) 

Consequently the Bergmann-Thompson gravitational energy-momentum pseudotensor 
density, found from a decomposition of the form 7bt = \g\^T + £bt, is 

Ct = V s - {glh^^p^ 6 . (30) 

In normal coordinates to zeroth order the pseudotensor £bt vanishes so Tbt = 
which is the desired limit inside of matter. In vacuum the complex ( J29i) reduces to the 
pseudotensor, which has its first nonvanishing contribution at the second order: 

2«c T = \g\ h [2T\ a v^ - 2r A /r^ A + 9 ^t X(7T t° Xt ] ~ ^|^(7£ + ^s)^- (3i) 

The better known Landau-Lifshitz expression (15] can be obtained from this 
transvected superpotential with an additional improper density weight factor ([7]): 

2kT£ := d x (\gfig" 5 Ug x 5 ) = 2k\ 9 \^ - |^/TV£%(|^^). (32) 
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A short calculation shows that this extra weight factor makes no contribution to the 
zeroth-order matter-interior limit nor to the second order result in vacuum. Thus the 
Landau-Lifshitz vacuum result is given by ( 13T1) . as was noted some time ago [23] . 

By the way, Goldberg [17] has proposed an infinite class of arbitrary density weight 
energy-momentum expressions obtained from those of Einstein and Landau-Lifshitz in 
a manner similar to that just considered. For essentially the same reasons, Goldberg's 
weighted density factors will not yield any modified values in the small region limits we 
are considering here. 



4-3. Papapetrou 

For our analysis of the Papapetrou [12] and Weinberg [18] expressions we use the Einstein 
tensor expansion 

7GT = 5Z& b l\9 ef g mn - 9 me 9 nf ] d aP g ef 

+ {2gP» g v - g^gP) (d a g a5 T s ^ - d^g aS T 5 ^ + T a Sa T% y - r a Sj r 6 Pa ) . (33) 

For the Papapetrou energy-momentum complex [TJl [T3], [H] we find, using (13"3l and the 
appropriate superpotential (JHJ) 



2„rr := 5Z& b d a p(\g\h mn ) (34) 

= W b (g ab - ^g ab )d a [\g\k\g ef g mn - g me g nf )d,ge f ] (35) 

~ 2\g\*GT + | 5 |i(2/V y - <T/ 7 ) [^r,^ + r W] 

+ 5zC^g ab \g\ h g me g nf d^ge f 

+ OSrtl* {-\g ei g f] g mn + g mi g e3 g nf + g me g m g f3 )d a g l3 d p g ef , (36) 

where we use for any metric expression AF := F — F. Here we have indicated one way to 
arrange things so that we can get the terms of the desired orders. In normal coordinates 
to zeroth order within matter the expansion (1361) clearly reduces to the desired result. 
In vacuum to second order a lengthy calculation yields 



XX 



2Kt» v -—[AB-S-K]^. (37) 
4-4- Weinberg 

With the help of (133|) and the appropriate superpotential (jSJ) for the Weinberg [T8] 
expression (note: the same expression is considered in [I] §20.2) we find 

2kt3T := 6zC\g\"g ab (-g me g nf + lg mn g ef )d a pg ef (38) 



d ma d nb\g\ 2 g [~9 9 J + ~9 9 )°a 



2 ^ a ,-«p9ef 



- 0#A 



\g\^g ab {-g me g nf + \g mn g ef ) 



d a/ 3g e f (39) 



2\ g \*G»» + \g\H^g lu - g^g^) + r a 57 rV] (40) 



Energy-momentum density in small regions: the classical pseudotensors 9 

- [aw + ,-vw + 9 -w - ±«-W')] ^ w 

Here we have indicated one way to get the desired orders. In normal coordinates, as 
expected, to zeroth order it gives the desired material limit; in vacuum to second order 
after a long calculation we found 

2Kt^^^[-B-2S-3K} t f\ (41) 

4.5. M0ller 

M0ller's holonomic energy-momentum complex [TU] follows from the superpotential (|10p . 
We find 

2kT» u := - <9 A (M V'TV)^ = d x \\g\h^9 X6 (d p g Sv - d s gp v )} (42) 



= \g\*g^g Xd (d xmv - d X5 g p „) + d x [\g\* g^ g xt> j {d p g Sv - d s g Pv ) (43) 
~ \g\*Rf> v - ^g sx {d X5 g Pv - d xp g 5v + d 5v g xp - d Pu g xs ) + 0(T 2 ) (44) 

~ \g\*RP v + 0(x 2 ). (45) 

Note that the M0ller energy-momentum complex fails to give the correct small region 
material limit (by the way, it also fails to agree with the proper asymptotic limit 
presented in p], §20.2). Using Einstein's equation and <$fy we have to zeroth order 

2t^ = -(T^ + ^T 5 s ), (46) 

which is non-vanishing inside of matter. Therefore the M0ller expression does not satisfy 
the equivalence principle. From our perspective Miller's holonomic expression is thus 
disqualified as a satisfactory description of energy-momentum. 

Although it is consequently only of academic interest, nevertheless, for 
completeness, we briefly report here on the small region vacuum limit of M0ller's 
pseudotensor. From a long calculation, which requires the 4th order Riemann normal 
coordinate expansion for the metric |30j . 



1 2 £ 

ga/3,a\^u = P( — — Raal3\; f iv + ^T-R a/3\R^au), (47) 

where P is the aXfiu symmetrization projection operator, we found to second order 

2Ktl v ^f^[2B- l -S-K] i f v . (48) 

5. Combinations of the classical pseudotensors 

The M0ller expression is ruled out of our further considerations by its material and 
asymptotic limit. The other four expressions are satisfactory asymptotically and in the 
material limit. However none of them give a good small vacuum value. Following the 
example of [23J we shall consider linear combinations of the four satisfactory classical 
pseudotensors. In the aforementioned work it was noted that in vacuum a certain 
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combination of the Einstein and the Landau-Lifshitz (or equivalently the Bergmann- 
Thompson) expressions is to second order proportional to the Bel-Robinson tensor; 
specifically they found 

B a /3fiu = d^ v (t E + |tLL)o/3- (49) 

This raises a couple of issues not addressed in that earlier work. In particular one 
should pay attention to the overall normalization in order to satisfy the small region 
matter interior limit (the same normalization will also give the correct magnitude for the 
total energy momentum of an asymptotically flat space). One also must cope with the 
mixed index positions which naturally occur for the Einstein pseudotensor in contrast 
with the two contravariant indices of the other three pseudotensors. Algebraically at 
the RNC origin it matters not whether one raises the index on the Einstein expression 
with g a P or g a " . However only the latter choice is allowed if we wish to preserve the 
conservation property. Thus we see another reason why we are paralyzed unless we 
allow for the explicit appearance of the Minkowski metric in RNC. Hence we define 

C := (50) 

With this convention we can consider the linear combinations (suppressing indicies 
for simplicity) 



2nt := 2n[et E + bt BT + pt P + wt w ] ~ (e + b + p + w)2G 



(51) 



1 

H — xx 
9 



7 

B(2e+ -b + Ap 



w) — S(—e — -b + p + 2w) 



K(p + 3w) 



where we use ~ to mean that only the zeroth order and second order vacuum values 
have been indicated. In order to have the correct asymptotic and material limits one 
must require e + b+p + w = 1, and for a good small vacuum limit to Bel- Robinson the 
coefficients of S and K should vanish — thus three restrictions on the four parameters. 
A convenient way to parameterize the set of acceptable coefficients is 

e = l-A, 6 = 2/3, p = 3A/2-l, w = -A/2 + 1/3, 



(52) 



which yields 



t(X) 



A)£e + g^BT 



A 1 . 

o - o ( 3 *p 



v-1 



+ A 



-t E + - (3tp — tw) 



(2k)" 1 [2G + XBxx/2] = T + XBxx/4k. 



(53) 

(54) 
(55) 

The combination in the first bracket of ( j54j) gives both the desired asymptotic and small 
region material result, but vanishes to 2nd order in the small vacuum region limit; the 
combination in the second bracket vanishes both asymptotically and in the material 
limit, while in the small vacuum limit it gives the desired 2nd order Bel- Robinson 
contribution. Formally we can choose any value for A. Physically we want A > 0. We 
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know of no principle to fix the magnitude of A. Some values stand out: A = 2/3 gives 
the (here properly normalized) case found earlier [23J: 




Note that the earlier work found just one expression. The restricted form of the 
expressions considered there (i.e., not explicitly containing g) are not in our opinion 
justified. To take linear combinations we need to get the indicies on the Einstein 
pseudotensor at the same level, that requires g^ u . Moreover, the Landau-Lifshitz 
expression is actually of the wrong density weight (this can be adjusted by including a 
suitable power of \g\). To the order considered in that work one can equally well use 
the Bergmann-Thompson expression, which we have used here. 

Comparing with the earlier work, we note that both the Weinberg and Papapetrou 
expressions cannot be constructed without explicitly using a reference. We certainly 
want to allow for these expressions — especially the latter which has been generalized to 
define the total energy for an asymptotically anti-de Sitter space [3l] . 

6. Conclusion 

In this work we found the values for the gravitational energy-momentum density given by 
the famous classical pseudotensors: Einstein, Papapetrou, Landau-Lifshitz, Bergmann- 
Thompson, Goldberg, M0ller, and Weinberg, in the small region limit to lowest non- 
vanishing order in normal coordinates. All except Miller's were found to have the zeroth 
order material limit required by the equivalence principle. For small vacuum regions 
we found that none of these classical holonomic pseudotensors satisfies the 2nd order 
vacuum criterion of being proportional to the Bel-Robinson tensor. 

However certain linear combinations of these classical holonomic pseudotensors do 
give the desired 2nd order vacuum result. Generalizing an earlier work [23] which had 
identified one case, we found another independent linear combination satisfying this 
requirement — and hence a one parameter set of linear combinations of the Einstein, 
Bergmann-Thompson, Papapetrou and Weinberg pseudotensors with this same desired 
property. (This had been overlooked because the earlier work had not allowed for the 
explicit use of the Minkowski metric as a reference in the pseudotensor.) 

Regarding the physical meaning of these combinations, although that has not been 
worked out in detail, there is a straightforward way to find it. Each viable parameter 
choice fixes a specific Hamiltonian boundary term — and thereby not just the value 
of the Hamiltonian but also an associated boundary condition. Nevertheless these 
combinations may still appear rather artificial, as being mathematically and physically 
contrived. This situation may be contrasted with the case of another (nonholonomic) 
energy-momentum expression: elsewhere we have shown that the tetrad-teleparallel 




t(2/3) = (l/3)t E + (2/3)t B T — T + Bxx/6k. 
The choice A = 1 gives another simple case: 

t(l) = (2/3)t B T + (1/6) (3t P — %) ~ T + Bxx/Ak. 



(56) 



(57) 
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energy-momentum gauge current expression [32] naturally has the desired Bel- Robinson 
property [33] . 

Although it still allows for a certain amount of freedom, the small vacuum region 
Bel-Robinson positivity requirement provides a strong restriction which excludes many 
otherwise satisfactory energy-momentum expressions. 
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